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Abstract 

The possibility to derive an equation for the mean velocity field in 
turbulent flow by using classical statistical mechanics is investigated. An 
application of projection operator technique available in the literature is 
used for this purpose. It is argued that the hydrodynamic velocity defined 
there, in situations where the fluid is turbulent, is to be interpreted as 
the mean velocity field; in that case, the momentum component of the 
generalized transport equation derived there is the mean velocity equation. 

In this paper, stationary incompressible flow for constant mass den- 
sity and temperature is considered. The stress tensor is obtained as a 
nonlinear functional of the mean velocity field, the linear part of which 
is the Stokes tensor. The formula contains a time correlation function 
in local equilibrium. Presently, there exists a microscopic theory for time 
correlations in total equilibrium only. For this reason and as a preliminary 
measure, the formula has been expanded into a power series in the mean 
velocity; though this limits the applicability to low Reynolds number flow. 
The second order term has been evaluated in a former paper of the au- 
thor. For the third order term, the form of the kernel function is derived. 
Its calculation with the aid of the mode-coupling theory is completed; it 
will be reported in an separate paper. An numerical application with the 
data of the circular jet is under way. 

Key words: statistical thermodynamics, projection operator technique, 
turbulent flow 

1 Introduction 

In turbulence experiments, very often the relevant variables - those which can 
be related to the external conditions of the experiment - are not the actual 
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quantities but their mean values. The process is then described by a statistical 
theory. It is well known that it can be very difficult to derive equations for mean 
values which are closed. It can be argued that the reason is that it is difficult 
to construct the multi-point probability distribution for the process. 

In this paper, the possibility to derive an equation for the mean velocity field 
in turbulent flow by using classical statistical mechanics is investigated. Then, 
in principle the probability distribution of the process can always be constructed 
starting from the total equilibrium distribution. On the other hand, it will be 
necessary to distinguish between macroscopic and microscopic parts of the 
motion and to formulate the latter in a way suitable for hydrodynamic purpose. 

It is known that the Navier-Stokes equation can be derived from Statistical 
Mechanics. This has been performed first by Chapman and Enskog in 1916 and 
1917, with the Boltzmann equation as a starting point; see, e. g., Huang [1]. 
More recently, a formalism has been developed which starts from the Liouvillc 
equation and applies projection operator technique; see Zwanzig [2], Mori [3]. 
Here, the presentation of Grabcrt [4] is used as a reference. In the relevant 
part of this work, a generalized transport equation is derived; the application 
for simple fluids and a suitable approximation of the stress tensor lead to the 
Navier-Stokes equation. In the present paper, arguments are given that for 
situations where the fluid flow is turbulent, the momentum component of the 
generalized transport equation for simple fluids actually is the mean velocity 
equation. The formulation is restricted to stationary flow in incompressible 
fluid of constant density and temperature. 

In section 2, some definitions of Statistical Mechanics and the basic steps 
which lead to the generalized transport equation are referenced from [4]. The 
interpretation of the momentum component as the mean velocity equation is 
explained. It is known that the form of this equation is very similar to the 
Navier-Stokes equation, with an additional friction force term (Reynolds force). 
There are numerous approaches to formulate this quantity (see, e. g., [5] and 
the references therein) ; it is seen that it is nonlinear in the velocity and, after all 
evidence, also nonlocal. - As a preliminary measure, the friction force has been 
expanded into a power series in the velocity; though this limits the application 
to low Reynolds number flow. This is explained in section 3. The main results 
for the 2 nd order term reported earlier [6], [7] are quoted; the calculation of the 
formula for the 3 rd order term is presented in this paper. 

2 Mean velocity equation 

In this section, the definitions and the basic steps of the derivation of the gen- 
eralized transport equation are referenced from Grabert [4]; the notation is 
slightly different. The fluid is considered to be a system of N particles of mass 
m with positions yj and velocities Vj (simple fluid) which are combined to the 
phase space vector z. Vector components are described by Latin indices, e. g. 
y rj = {yja} ■ The particles interact with a symmetric, short-ranged, pairwise 
additive interparticle potential. The system is enclosed in a box of Volume V . 
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A function g(z) is called a phase space function, or microscopic variable. Espe- 
cially, we need the space densities of the conserved quantities particle number, 
energy and momentum n, e, p which are collected to a 5-element linear matrix 
a . They are functions of an additional space variable x: 

N 

a = ^a j 5{x-y j ) (2.1) 

3=1 

For the particle functions a d - we have rij = 1 , Pj = mVj , while the en- 
ergy function contains the interparticle potential. The quantities a obey the 
conservation relations: 

a = -Vs (2.2) 

The fluxes s have the same general structure as the a (2.1); the particle 
functions can be found in [4]; especially, we have Si = p/m . The time evolution 
of any phase space function g is described by the Liouville equation: 

9 = i£g (2.3) 

iC is the Liouville operator, a linear differential operator the form of which 
can be found in any textbook of statistical mechanics. From (2.3), the formal 
solution for g(t) given the initial value g is: 

g(t)=e iCt g (2.4) 

In the statistical model, z and N are considered random variables; that is, 
the probability density f{z, N) is of grand canonical type. The ensemble mean 
value (expectation) of a phase space function g is defined in the 'Heisenberg' 
picture: 

(»)(*)= E / dzg(z,N,t)f(z,N) (2.5) 

N=l 

In this formula, f(z, N) is the initial probability distribution, and g(z, N, t) 
is the value of g at time t if the initial positions and velocities of the particles 
are described by z . The operation (integration + Summation) is sometimes 
indicated by the symbol 'tr': 

oo „ 

tr{il} = E / d*to(*> N ) ( 2 - 6 ) 

jv=i J 

Certain probability densities (also called distributions here) are frequently 
used in the analysis. One of them is the (total) equilibrium distribution which 
corresponds to macroscopic rest: 

fo - exp(<f> + f3(nN - H{z))) (2.7a) 

m = jpSj fN (2Jb) 
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Here, h is Planck's constant, — l/(ksT) , fcs being Boltzmann's constant 
and T the temperature, \i is the chemical potential which is a function of mass 
density and temperature, and H{z) is Hamilton's function which describes the 
total energy of the fluid. For the normalization constant, we have <J>o = —/3PV, 
P being the equilibrium pressure. Expectations with respect to the equilibrium 
distribution are denoted by ()o- In case of a simple fluid, the 'relevant probability 
distribution' of Grabert's formalism (see [4], sec. 2.2) is the local equilibrium 
distribution: 

f L (t) = i>(N) exp($(t) - a(z) * b(t)), (2.8a) 
b = {/3(yu 2 - M ),/?,-/3«}, (2.8b) 
$(i) = -log(tr{V>exp(-a*6(i))}). (2.8c) 

Here the symbol * is introduced for the operation: Multiplication, plus Sum- 
mation over the 5 elements of the linear matrices a , 6, plus Integration over 
geometrical space. The elements of b arc called the conjugate parameters; they 
are functions of the quantities (3 , /i and u which we will sometimes call the 
thermodynamic parameters. [3 = l/{kT), with k being Boltzmann's constant 
and T the absolute temperature; fi is the chemical potential which is a function 
of temperature and pressure, and u is the hydrodynamic velocity. These quan- 
tities will be considered slowly varying functions of space and time. The b are 
defined such that the expectations of the a are identical to their expectations 
in local equilibrium: 

(a) = (a) L (2.9) 

The projection operator techniqe (POT) is a means for separating macro- 
scopic and microscopic parts of the random variables. It starts by defining the 
set of phase space functions which are relevant for the description of the process. 
For simple fluids, this set is identified with the densities of conserved variables, 
a . A projection operator is defined which projects out of any microscopic 
variable g the part which is proportional to the relevant variables. It reads: 

Vg = (g) L + (g5a) L * (8a <5a>^ * Sa (2.10) 

Here, Sa = a — (<i}l; ()2 denotes the inverse of the expectation matrix in 
the formula. For stationary flow, V is time independent. The analysis in [4] 
starts by splitting the exponential operator (2.4): 

e lCt = e lCt P+ f dt 1 e lCt ' PiLil-P)^ 1 -^ 1 ^-*'^ +{1-P)c {1 - V ^ ct (2.11) 
Jo 

This corresponds to [4] , formula (2.4.1), specialized to stationary flow, where 
especially P =const(i). With (2.11), the Liouville equation (2.3) is reformulated: 

a = c tCt PiCa+ 

+ / dt'e tCt ' PiC(l-P)c^ v ^ t - t 'UCa+{l-P)c (1 - v ^ ct iCa (2.12) 
Jo 
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By averaging over the initial probability density, and after some manipula- 
tions, Grabert's generalized transport equation [4], (2.5.17) ist obtained. Below, 
this equation is presented for stationary flow in simple fluids. It is postulated in 
POT that the initial probability density is of the form of the relevant probability 
density. Grabert states that this should not be considered a general restriction 
of the method but a means to form the general particle system into the type 
specially considered; see [4], sec. 2.2 . In [4], sec. 8.3, it is shown that for sim- 
ple fluids the relevant probability density is that of local equilibrium. For the 
present approach this means that turbulent flow is considered which initially de- 
veloped from laminar flow with suitable velocity gradient. - It is a consequence 
of this postulate that the last term in (2.12) vanishes after averaging. Moreover, 
it is shown that: 

{e tCt ViCa) = (a J = -V • (s) L (2.13) 

In the last step, the conservation relations (2.2) are introduced. Station- 
ary flow is considered to be the process described by the generalized transport 
equation under stationary conditions and for very large times. One obtains: 

= -V-(s) L + D (2.14) 

D a (x) = -W c J dx' J ijleP-^'^.^H^^lJviWa;') (2.15) 

s ac {x) = {1 - V)s ac {x) (2.16) 

Latin and greek indices run over 3 and 5 values, respectively; this is some- 
times expressed by saying that, e. g., the index a 'runs over 1, 2 and the latin 
index a\ Eqs. (2.14), (2.15) correspond to [4], (8.1.13), (8.5.1) for stationary 
flow. From these, now the momentum equation is taken, for which one obtains 
([4], (8.4.15)): 

(Sac) L = PU a Uc + P&ac (2.17) 

Here, S ac is the Kronecker symbol, p the mass density and u the fluid velocity 
defined by: 

(p)=pu (2.18) 

([4], (8.3.12)). At this point, the continuity equation is introduced. This ac- 
tually is the mass density component of (2.14), and for incompressible constant 
density flow it reduces to: 

V-m = (2.19) 
With these formulas and (2.8b), (2.14), (2.15) obtain their final form: 

pu ■ Vit = -VP + V • R (2.20) 

Rac(x) = J dx' ')V d u b (x') (2.21) 

S abcd (x,x')=p dt{[e^- v ^ Ct s ac {x)]s bd {x') L (2.22) 
Jo 
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It is seen immediately that (2.20) is the hydrodynamic velocity equation for 
stationary incompressible flow. The stress Tensor R is, in general, a nonlinear 
functional of u. Grabert, in [4], sec. 8.5, performs an approximation to the first 
order of Vu and obtains exactly the Stokes form of the tensor. Therefore, in 
this approximation, (2.20) is the stationary Navier-Stokes equation. 

By (2.18), u is essentially equal to the expectation of the momentum density. 
By definition from probability theory, the expectation ist an average built from 
a set of realisations of the process, which can be interpreted as repetitions of the 
experiment und identical external conditions. If the flow exhibits macroscopic, 
i. e. turbulent, fluctuations, the averaging process includes these. There is 
barely another possibility than, in this case, to interpret u as the mean velocity 
of the flow. On the other hand, if one wants to define the 'point' velocity of 
turbulent flow, it would be necessary to introduce a conditional expectation, 
which excludes the macroscopic fluctuations, with respect of which the mean 
quantity would still be a random variable. - The preceding statement is quite 
general; it is still valid if one performs a projection operator analyses. Thus, in 
case of stationary turbulent flow (2.20) is the mean velocity equation. It is of 
course necessary to check this theoretical statement by bringing the equation 
into a form that can be evaluated, and comparing the results with a turbulent 
flow experiment. 

The difference is that for turbulent flow the nonlinear part of the stress 
tensor R is essential. The definition formula for the stress tensor kernel function 
S (2.22) contains a time correlation function in local equilibrium. This is a 
quantity which, for processes with constant mass density and temperature, is a 
functional of the velocity field. It has to be evaluated in advance, by a separate 
statistical-mechanical formalism. At present there is no theoretical means to 
perform this; instead, it is possible to calculate correlation functions for total 
equilibrium. It should be emphasized that, since total equilibrium corresponds 
to macroscopic rest, the latter quantities do not depend on the flow properties; 
they are material 'constants' of the fluid. - In order to pursue the analysis, S has 
been expanded into a functional power series in u; as will be seen, the coefficients 
of the series contain total equilibrium correlations. It will be necessary then to 
evaluate the lowest order terms of the expansion and to work with formulas 
(2.20), (2.21) with an approximated quantity S. When the formlism is applied 
to a given flow configuration, it is possible to render all variable quantities in the 
expansion dimension-free, which causes certain constant factors to appear in the 
terms of it. The author elaborates on an application to circular jet flow; where 
these factors show up actually as increasing powers of the Reynolds number Re. 
Thus, it is seen that in this case the application should be restricted to low 
Reynolds number flow, just beyond the laminar-turbulent transition. 

3 Expansion of the kernel function 
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S depends on u via the quantities b (2.8b) in the formula for the local equilib- 
rium distribution. The series expansion of S is performed by expanding it with 
respect to the b at the point b — 6 which corresponds to u =0: 

b = {-/3 M ,/?,0} (3.1) 

b-b o = {^u\0, -Pu} (3.2) 

Replacing b by 6o changes the local equilibrium distribution (2.8a) into the 
total equilibrium distribution (2.7a) with the prescribed mass density and tem- 
perature. The power expansion of S reads: 

S = S\b + S i\b a * (6 " bo) + |||| bo * *{(b - bo), (b - bo)} + ■■■ 

= S (0) + S (i) + S (2) + . . . (3 3) 

In the second row, symbols are applied to the different orders of the expan- 
sion. The expansion in u is obtained from (3.3) by restricting the summation 
inherent in the ^-operation to the last element of b bo hi (3.2). It is seen there 
that the first element also depends on u. But for the present investigation, 
where no sound or heat conduction processes are considered, the expansion co- 
efficients come out to be non-zero only for 'Latin' values of the indices; thus the 
first element of b does not contribute. 

The linear part of R which results from the constant term of the ex- 
pansion (3.3) leads to the Stokes form of the stress tensor, which coincides with 
the result of Grabert. At present, the 2 nd and 3 rd order terms of R have been 
analyzed. From (2.21) it can be seen that these stem from the parts and 
of the S-series, respectively. To obtain these formulas, it is necessary to 
calculate the first and second order functional derivatives of S. The first of 
these has been done in Piest [6]. For completeness, the derivation is repeated 
in the appendix. The result is formula (a. 17); by setting b = b one obtains: 

5S abcd {x,x') 
Sb e (x") lb=b « 

f-OC 

= -(3 / dt([e^-^ iCt (l - V )s ac (x)][(l - V )s bd (x')}p e (x")) (3.4) 
Jo 

Vo is the total equilibrium projection operator: 

'Pag = (g)o + (gS a) * (SoaSoa)^ 1 * S Q a (3.5) 

Here, Sog = g — (s) : an d we have (p) Q = 0.- The total equilibrium triple 
correlation function contained in (3.4) can be calculated using the mode coupling 
technique in the form of Martin et al. [8], which is a systematic version of 
the method of Kawasaki [9] , and has been complemented by Deker and Haake 
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[10]. The formulas have been evaluated by the author, and applied using the 
experimental data of the circular jet (Piest [7]). One obtains the solution again 
in form of a power series. The expansion parameter \ reads: 



1 



(3.6) 



X = 



px 



v is the kinematic viscosity, x the distance of the observation point from 
the orifice. With a typical value x = 0.3 [m], for a laboratory experiment in 
air, one obtains x — 5-8 • 10~ n , which is so small that practically only the zero 
order term of the expansion counts.- It was rather surprising that for the linear 
term of the it-series of S (3.3), a zeroth order %-Term is found which leads to a 
second order term of the friction force D = V • R: 



A = a/(-fCv) is a physical parameter of the fluid, a being the thermal ex- 
pansion coefficient, 7 the isothermal compressibility, and cy the specific heat for 
constant volume. The apparent existence of this term poses a problem to the 
intcpretation of equation (2.20). It should be noticed that, since for sufficient 
low Reynolds number the flow is laminar, there are no macroscopic fluctuations 
then, and u equals the classical hydrodynamic 'point' velocity. Therefore, for 
decreasing Re, equation (2.20) should reduce to the stationary Navier-Stokes 
equation. On the other hand, Z)*- 2,0 ' is of equal power in Re as the quadratic 
term on the left of (2.20) so that, even for small Re, (2.20) remains different. 
Since this is not possible, something must be wrong with the prerequisites of 
the derivation of (3.7). As a preliminary measure, it is assumed that the term 
actually does not exist. It is of course one of the most urging requirements for 
this approach to sufficiently explain this defect. 

The third order term of the expansion of R contains, after (2.21) and (3.3), 
the second-order derivative, which is presented in (a.35). When this formula is 
applied to b = 6 as in (3.4), the summation in the corresponding term in (3.3) 
is restricted to the last term of (3.2) and the result is inserted into (2.21), on 
obtains for the friction force D = V • R: 



(x) = / dx'dx"dx"'K abcd (x, x', x", x'")u b {x')u c {x")u d {x'") (3.8) 



For formal reasons, all V-operations have been transferred to the correlation 
function, partly by partial integration. Formulas (3.8), (3.9) are the main result 
of the present paper. The definition formula for the kernel function K contains 




(3.7) 



1 r°° 

K abcd {x,x',x",x'") = --/3 3 / dt([c (1 - V ^ Ct V e (l -V )s ae (x)}x 



x [V f (l - V )s bf (x')}5o(p c (x")p d (x"'))) 



(3.9) 
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a time integral over a correlation function which is double in time but quadruple 
in space. To calculate it, the author has again applied the technique described 
in Martin et al. [8], Deker and Haake [10]; it had to be enlarged slightly so that 
quadruple correlations could be determined. This investigation will be presented 
in a separate paper. The resulting formula is rather lenghty and will not be given 
here. It should be emphasized that with K obtained as an explicit formula of 
the four space variables, (2.20), (3.8) form a closed system for calculating the 
velocity field for a given flow configuration. Moreover, since this formula is the 
lowest order term of the expansion (see the comment to (3.7)), (2.20), (3.8) is 
the simplest form of the system for checking, by comparing with experimental 
results, whether the approach works. As has been mentioned, a numerical test 
using the data for the circular jet is in progress. 

4 Summary 

An approach to arrive at the mean velocity equation for turbulent fluid flow has 
been attempted with the aid of the projection operator technique in classical 
Statistical Mechanics. The hydrodynamic velocity is defined in this technique 
via the conjugate thermodynamic fields in the formula for the relevant proba- 
bility density; multiplied by the mass density, it is identical to the expectation 
of the microscopic momentum density. It is argued that in situations where the 
fluid flow is turbulent, this is precisely the mean velocity field of the flow. If 
this argument is correct, the momentum component of the generalized transport 
equation derived by this technique is the mean velocity equation. 

Stationary incompressible flow for constant mass density and temperature 
is considered. The formula for the stress tensor is a nonlinear functional of 
the velocity, the linear part of which has the form of the Stokes tensor. The 
formula containes a local equilibrium time correlation function. At present, 
there exists a theory for calculating correlation functions for total equilibrium 
only. As a preliminary measure, the stress tensor has been developed into a 
power series in the velocity, though this limits the applicability of the equation 
to low Reynolds number flow. The coefficients of the expansion contain total 
equilibrium correlations which can be calculated. The second order term has 
been evaluated in a former paper of the author. For completeness, the main 
results have been reported here. The calculation leads to a second order friction 
term which is comparable in form to the convolution term of the equation. This 
constitutes a problem to the present approach since from general knowledge 
about the Reynolds equation such a term cannot appear. 

For the third order term, the form of the kernel function has been calculated. 
The formula has been evaluated with the aid of the mode coupling theory of 
Statistical Mechanics; the results will be reported in a separate paper. A nu- 
merical calculation in order to test the equation with the experimental data for 
the circular jet is under way. 
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Appendix: 

Calculation of functional derivatives 

In this appendix, the first two derivatives, with respect to b, of the kernel 
function S (2.22) are calculated. The right hand side of (2.22) depends on b at 
four different places; thus we may write: 

SS abcd (x,x') _^5S® 

6b e (x") ^ Sb 1 j 

The definition of the four terms is given in the formulas to follow. In addition, 
certain auxiliary formulas are written which can be verified directly. For the 
first term, we use the rule: 

^ = -hSa (a.2) 

The first term reads: 



5S_ (1) a [°° JMM __ ( Sf, 
~8b 



3 1 dt tr{ Jb^) [e(1 ~ V)tCt '"W*" 



= -f3 dt([e^ lCt s ac (x)]s bd (x')6p e (x")) L (a.3) 
Jo 

Next, we need the formula for the differentiation of the projection operator 

V : 

¥L = -V5a(l-V) (a.4) 

The second term reads: 



r°° A-p 
-15 jf dt{[eV-rW ls ^8 ac (x)]8 bd {x')) L 



5S& 

Sb 

r OO 

= dt([e^ iCt Vs ac (x)5p e (x")}s bd (x')) L (a.5) 
Jo 

When (a.4) is applied to the third term, one finds, after some manipulations: 

f (3, = - /? ^ <it ( [ ea-^«^ wl _g_ SMM)l 

= d \ dt([(l - V) s ac (x)]Vs hd (x')5 Pe (x")) L 
Jo 

= (a.6) 

In the second step, it is used that we have: 

e {l-V)iCt^ _ -p) = (! _ -p) e (l-V)i£t^ _ ( a 7 ) 
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([(1 - V) 9l ]Vg 2 )L = (a.8) 

(a. 8) is valid for any microscopic functions g\, g 2 ■ - The part containing the 
differentiation of the exponential operator is a somewhat more involved. The 
differentiation formula is: 

i P (i-p)ia rt , a 

/ ,./ CI _T>'W/V — ,. a 



Sb 



[ dt'e^-^t'pSa^e^Wt-t') (aQ) 
Jo dt 



For the corresponding time integral, it follows after interchanging the suc- 
cession of integrations: 

/ dt dC = / dt' dteV-TWvSpeix'^eQ-W* 

J 6b f{x'») J Jo F ' 'dt 

= / dt'{\im e^ iCt ' VSp e ix")e^ iCt -e^ iCt ' V6p e ix")} (a.10) 

Jo t^oo 

One obtains for the fourth term: 

=/3 i dt{[ 5b e (x») *™Wu&)h 



roc 

/ dt'{ lim ({e^-W V S Pe (x") e ^ lCt s ac (x)]s bd (x')) 

Jo t^oo 



L 



([e (i- W v 6 Pe (x")s ac (x)]s bd (x')) L } (a.ll) 



We have to evaluate the limit expression. We will assume here that in the 
limit of large times, the factors of a time correlation will become statistically 
independent so that for stationary processes we have: 

lim {[e^ iCt A]B) L = (A) L (B) L (a.12) 
For any phase space function g, let us define a quantity F(t'): 

F[t') = lim ([cd- W V 5p e ix") e^-W s ac ix)]g) L (a.13) 

t— >oc 

The evaluation results in: 
Fit') = \^{{5 Pe ix")^-^ t s ac ix)) L {g) L 

t — >oo 

+ (6p e ix")[e^ iCt s ac ix)]6a) L * (SaSatf * ([e^ iCt ' Sa]g) L } 
= (s ac ix)) L {(5 Pe ix")) L (g) L + (5p e ix")5a) L * 

* (SaSa)^ 1 * ([e^-W 5a]g) L } = (a.14) 

The first step is the evaluation of V using (2.10); next, we use (a. 12); finally, 
the factor outside the curled brackets vanishes, since we have, for any phase 
space function g: 

<(1 - P)S>l = (a.15) 
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Therefore, in (a. 11), the limit term vanishes, and we have: 

-£ =P J dtfeV-TW Vs ac (x)Sp e (x")}s bd (x')6p e (x")) L = 

(a.16) 

Thus, we obtain: 

SS abcd (x,x') = SS_ W 
Sb e (x") Sb 



= -0 dt([e^Ws ac (x)]s bd (x')Sp e (x")) L (a.17) 
Jo 



The calculation of the coefficients of the second order derivative of S parallels 
to a certain extent that of the first order. The right hand side of (a. 17) depends 
on 6 at five different places; thus we write: 

S 2 S a bcd(x,x') _ \U S 2 S ^ 



5b e (x")6b f (x'") f^iSbSb 



For the first part, we use (a. 2): 



(5 2 S (1) 



5b Sb 



= -/3^dttr{^^[e( 1 -^ t ? ac ( : r)]^( a; ') ( 5p e ( a; '')} 

= (3 dt([e^ iCt s ac (x)]% d (x')5 Pe (x'')6p f (x''')) L (a.19) 
Jo 



The second term , with (a. 4), turns out to be: 



r2 o (2) /*OC r-p 

0J ^([e( 1 -^)^ ^^ Sae ( a; )]g bd (x0^e(x ,/ ))i 



= -/?/ ^([e( 1 -^ t 7's ac ( a; )%( a; '")]s M ( : r')^e(a ; "))L (a.20) 
Jo 



'0 

In the same way, we have for the third term: 



r2 o (^) />oc r/p 

' = /?j( d t([e( 1 -^J Qc(a; )]^-^ SM ( a; ')^ e ( a; ''))L 

/•OO 

/?/ dt([e( 1 - T ')^ t s ac ( : r)][Ps M ( a; ')%(^")]^e( : r"))L (a.21) 
Jo 



56(56 



For the fourth part containing the differentiation of the exponential operator 
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we have, with (a. 10): 

S 2 S (4) f°° s e (i--P)iCt 

SbSb = - P J dt{[ Sb f (x»') ^{x)]s bd {x')8p e {x")) L 

= -(3 / ^'{lini([e( 1 - T ')^>%( a 3"')c( 1 - T, ) 4£t 5 ac ( : r)]5 M ( ; E')^ e (^"))L 
Jo 

- ([c( 1 -^ ct >%( : r / ")s ac ( a; )]s bd ( a; > 13e ( a; "))L} (a.22) 

The evaluation parallels that of the corresponding term of the first derivative. 
We finally find: 

c2 ^ ( 4 ) f-OO r2 o ( 2 ) 

(3 / dt([c( 1 - T ')^7's ac ( : r)%-( a; "')]s bd ( a; ')^ e ( a; "))L ' 



6b8b ^ J u ' /J"~v~ //- Jb56 

(a.23) 

For the fifth term, we need the differentiation rule for (a )l- 

5 M± = -(SaSa) L (a.24) 
oo 

We obtain: 

P OC 

= -/?/ rfi([e( 1 - p ) i£t ? oc ( a; )]? 6d (a ; '))L<^e(a ! /, )%(a ; / ")>L (a-25) 
Jo 

Finally, we want to show that the term (a. 21) vanishes. The first step is to 
evaluate V with the aid of (2.10): 



6>S (3) 



= -(3 / cit{(s M ( : r')^ / ( a; "')}L([e (1 - T ' )4£ *s QC ( : r)] ( 5p e ( : r"))L 
Jo 



Sb 5b 

+ (s bd (x')Sp f (x"')Sa) L * (SaSa)^ * {[e^-^ iCt s ac (x)]SaS Pe (x")) L } (a.26) 

The first term on the right vanishes. Both factors are zero; e. g.: 

(s bd (x')Sp f (x'")) L = (s bd (x')(l - V)Spf(x'")) L = (a.27) 

To investigate the second term of (a.26), we introduce some auxiliary func- 
tions: 

Z= / dte^- r ^ ct s ac (a.28) 
Jo 

We have, with a suitable chosen Y: 

Z = (1 - V)Y (a.29) 
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Therefore: 

(Z) L = (a.30) 

and: 

(Z5a) L ={) (a.31) 

Moreover, we write: 

S = (ZSaS Pe ) L (a.32) 
We want to show 5 = 0. We consider the identity: 

S(ZSa) L 5Z 5(a) L 

-^b^ = { sb- e Sa)L - {z)L ^b7 + * (a ' 33) 

The left hand side is zero because of (a. 31); so is the second term on the 
right because of (a. 30). Moreover: 



= r dt ^ e{i ~ v)ict ^- e(i ~ v)ict Z s ^ 

0b e J Obe Obe 



f 

Jo 



dt 1 lim c^- v ^ ct ' V 5p e e^- v ^ ct s ac {x) (a.34) 

t— >oo 



The first step is by direct calculation. To the terms, we apply (a. 9), (a. 4) 
respectively, to obtain the second step. Comparing this with (a. 13), we find that 
the first term on the right of (a. 33) is F(t') applied to g — 5a, and is therefore 
zero. Thus, we have shown that actually (a. 26) vanishes. In total, we obtain 
from (a.18), together with (a.19), (a.23), (a.24): 



2 o X 2 c, (5) 



S^S abcd {x,x') 5 Z S" S Z S 



5b e (x")6b f (x'") 5b 5b 5b 5b 

= dt{[e^ iCt s ac {x))s bd {x')x 
Jo 

x {5p e (x")5 Pf (x>") - (5 Pe (x")5p f (x'")) L }) L (a.35) 

References 

[1] K. Huang (1963) Statistical mechanics. Wiley, New York 

[2] R. Zwanzig (1960) Ensemble method in the theory of irreversibility. J 
Chem. Phys. 33: 1338 

[3] H. Mori (1965) Transport, collective motion, and Brownian motion. Progr. 
Theor. Phys. 33: 423 

[4] H. Grabert (1982) Projection operator techniques in nonequilibrium statis- 
tical mechanics. Springer, Berlin, Heidelberg, New York 



14 



[5] W. D. McComb (1994) The physics of fluid turbulence. Clarendon, Oxford 

[6] J. Piest (1989) Theory of turbulent shear flow. I. Kinetic theory derivation 
of the Reynolds equation: Avoiding the closure problem. Physica A 157: 
688 

[7] J. Piest (1990) Theory of turbulent shear flow. II. Calculation of the coef- 
ficient function. Physica A 168: 966; (1992) Theory of turbulent shear 
flow. III. Numerical test. Physica A 187: 172 

[8] P. C. Martin, E. D. Siggia, H. A. Rose (1973) Statistical dynamics of 
classical systems. Phys. Rev. A8: 423 

[9] K. Kawasaki (1970) Kinetic equations and time correlations functions of 
critical fluctuations. Ann. Phys. 61: 1 

[10] U. Deker, F. Haake (1975) Fluctuation-dissipation theorems for classical 
processes. Phys. Rev. All: 2043 



15 



